Quantum interference from sums over closed paths for electrons 
on a three-dimensional lattice in a magnetic field: 
total energy, magnetic moment, and orbital susceptibility 

Yeong-Lieh LinB and Franco Nori0 
, Department of Physics, The University of Michigan, Ann Arbor, Michigan 48109-1120 

We study quantum interference effects due to electron motion on a three-dimensional cubic lattice 
in a continuously-tunable magnetic field B of arbitrary orientation and magnitude. These effects arise 
from the interference between magnetic phase factors associated with different electron closed paths. 
The sums of these phase factors, called lattice path-integrals, are "many-loop" generalizations of 
j /^ ' the standard "one- loop" Aharonov-Bohm-type argument, where the electron wave function picks 
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up a phase factor e 1 * each time it travels around a closed loop enclosing a net flux <!>. Our lattice 
path integral calculation enables us to obtain various important physical quantities through several 
different methods. The spirit of our approach follows Feynman's programme: to derive physical 
quantities in terms of "sums over paths" . From these lattice path- integrals we compute analytically, 
for several lengths of the electron path, the half-filled Fermi-sea ground-state energy, Et(B), of 
noninteracting spinless electrons in a cubic lattice. Our expressions for Et are valid for any strength 
of the applied magnetic field in any direction. Moreover, we provide an explicit derivation for the 
00 ' absolute minimum energy of the flux state. For various field orientations, we also study the quantum 

interference patterns and Et(B) by exactly summing over ~ 10 29 closed paths in a cubic lattice, 
each one with its corresponding magnetic phase factor representing the net flux enclosed by each 
path. Furthermore, an expression for the total kinetic energy Et(B,v) for any electron filling v 
close to one half is obtained. We also study in detail two experimentally important quantities: 
the magnetic moment M(B) and orbital susceptibility x(B) at half-filling, as well as the zero-field 
susceptibility x(aO as a function of the Fermi energy fi. 
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The quantum behavior of noninteracting tight-binding electrons on a two-dimensional (2D) lattice immersed in a 
perpendicular magnetic field has attracted much attention due to its important role in diverse areas of physics. For 
instance, the problem is intimately related to the quantum .Hall effect. During the past few years, the ground-state 
energy of this system has also been intensively investigatedcTu in connection with mean-field studies of the t-J and 
| Hubbard models of high-T c superconductors. At half- filling, the absolute minimum of the kinetic energy has a half 
t-H ■ flux quantum ($o/2) per plaquette — called the flux state. In addition, there are local minima at flux values equal to 
' $o/2"* with m > 1. 

More recently, several groups have paid attention to the three-dimensional (3D) case: the kinetic energy of a 3D 
noninteracting electron gas under the influence of both a strong periodic potential and a magnetic field. For instance, 
for several rational values of the magnetic field, Skudlarski and Vignaleu analyzed the changes in the ground-state 
properties induced by the addition of hopping in the direction parallel to the uniform magnetic field. HasegawaO 
studied the density of states and the total energy with the external field iftthc (0, 0, 1), (0, 1, 1), and (1, 1, 1) directions 
■ for certain selected rational values of the resulting flux. Kunszt and Zeefl showed that for rational values of the flux 
the eigenvalue problem can be reduced to a one-dimensional hopping in momentum space. They then calculated jthp. 
energy spectra and the density of states for various selected rational values of the flux states. These three worksoEl 
' focused on a discrete set of (rational) magnetic field values. Here we are interested in obtaining results that are 
^ , continuous functions of the applied field — valid for both commensurate and incommensurate flux values. 

In this work, we focus on the isotropic 3D cubic lattice, namely, the hopping integrals in the x-, y-, and z-direction 
are the same and taken to be equal to one. The system is described by the Hamiltonian 

H = y^cl Cj evp{iA tj ), (1) 
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where (ij) refers to nearest-neighbor sites, and the phase = 27r J? A-dl is 2ir times the line integral of the vector 
potential along the bond from i to j. Throughout this paper, the flux quantum $o = hc/e is set equal to one. The 
goal of this paper is to explore conceptually different viewpoints and approaches to study this interesting problem. 
Focusing on the hopping motion of electrons on the lattice, we first study the quantum interference between the phase 
factors of electron closed paths. This phenomenon is the source for the lowering of the total energy in a magnetic 
field. 
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In Sec. II, we examine the quantum interference effects originating from lattice path-integrals — defined as sums over 
magnetic phase factors on different electron closed paths. There, the physical meaning of these lattice path-integrals 
S21 and the techniques we use to compute them are discussed in detail. Results for S21 in an arbitrarily oriented field 
of any strength are also presented. Furthermore, we obtain the lattice path-integrals for several flux orientations by 
exactly summing an enormous number (~ 10 29 ) of closed paths, each one weighted by its corresponding acquired 
phase factor. 

In Sec. Ill, we present an analytical calculation of the total kinetic energy Et(B) of the half-filled Fermi sea oi 
tight-binding electrons. We use two different approaches: one. based on a direct series expansion for the energyo 
and the other on a moment expansion for the density of states.lD These two different methods yield the same results. 
Moreover, they provide different insights into the problem. For an arbitrary field orientation, analytic results for the 
total energy, in terms of the lattice path- integrals, are obtained. It is found that the lowest energy state is reached 
when the fluxes per plaquette on each one of the xy-, yz- and zs-planes are equal to &o/2. We also show the variations 
of the total energy Et as a function of the flux, for various orientations of the magnetic field. Furthermore, we obtain 
an analytic expression for the total kinetic energy Et(B, v) for any filling v close to one half. 

In Sec. IV, we investigate the magnetic moment M(B) and orbital susceptibility of this quantum system. 

For various flux orientations, we also obtain the zero-field susceptibility x(aO as a function of the Fermi energy /x. 
We show that the magnetic response, in the presence of a strong periodic potential, is significantly different from the 
familiar Landau diamagnetism. The zero-field susceptibility x(p-) as a function of the Fermi energy fi exhibits a large 
diamagnetic (i.e., negative x) response at very low electron filling. For increasing fi, increases and fluctuates 

around zero. Paramagnetism prevails for large /i. Indeed, the orbital response is paramagnetic at and near half-filling 
(// = 0). We obtain the field dependence of M(B) — which can be regarded as a generalized current in a multiply- 
connected geometry — and x(B) at half-filling, and both show oscillations between para- and diamagnetic behaviors 
as a function of the flux. The frequencies of these oscillations, as functions of the flux, tend to decrease for increasing 
field. 

In Sec. V, we summarize our results. 

II. QUANTUM INTERFERENCE FROM SUMS OVER CLOSED PATHS: LATTICE PATH-INTEGRALS 

A. Physical interpretation 

The lattice path-integral of order 21 is defined as 

S21 = £ e** r , (2) 

All closed 21— step 
lattice paths T 

where <&r is the sum over phases of the bonds on the path T of 21 steps starting and ending at the same site. Let \ipi) 
denote a localized one-site electron state centered at site i. It is not difficult to notice that £2; corresponds precisely 
to the quantum mechanical expectation value (i(!i\H 2l \ipi), which summarizes the contribution to the electron kinetic 
energy of all closed paths of 2/-steps. 

The physical meaning of S21 (= {i'iW 21 ^^) thus becomes clear. The Hamiltonian H is applied 21 times to the 
initial state \ipi), resulting in the new state H 2l \i(ji} located at the end of the path traversing 21 lattice bonds. Because 
of the presence of a magnetic field, a magnetic phase factor e lAij is acquired by an electron when hopping through 
two adjacent sites i and j. S21 is nonzero only when the path ends at the starting site. In other words, £2; is the sum 
of the contributions from all closed paths of 21 steps starting and ending at the same site, each one weighted by its 
corresponding phase factor e** r where <J>r/27r is the net flux enclosed by the closed path T. 

It is important to stress that <i>r depends crucially on the traveling route of the path. For instance, <i>r will 
be positive (negative) by traversing a polygon loop counterclockwise (clockwise). Therefore, quantum interference 
information contained in £2/ arises because the phase factors of different closed paths, including those from all kinds 
of distinct loops and separate contributions from the same loop, interfere with each other. Sometimes, the phases 
corresponding to subloops of a main path cancel. 

B. Analytical computational formalism 

We will now compute the lattice path- integrals £2/. This is a difficult task since £2/ involves an enormous number 
of different paths (growing rapidly when the order increases), each one given by its corresponding net magnetic phase 
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factor. We have considerably simplified this calculation by successively iterating the recursion relation and analyzing 
the symmetries of the problem. 

We consider a unit spacing for the cubic lattice. The vector potential of a general magnetic field B = (B x , B y ,B z ) 
can be written as 

A = l;(zB y - yB z , xB z - zB x , yB x - xB y ). 

Let a/27r, b/2ir, and c/2ir represent the three fluxes through the respective elementary plaquettes on the yz-, zx- and 
xy-planes. Thus, a flux configuration is specified by (a, b, c). 

From the definition of S 2 i, it is clear that: (i) S obviously equals 1, (ii) S 2 i+i are always zero because there is no 
path with an odd number of steps for returning an electron to its initial site, (iii) the <S2z's are gauge invariant, and 
(iv) in the absence of the magnetic field, £22 are just the total number of 21-step paths on the cubic lattice starting 
and ending at the same site. 

First, it is instructive to evaluate the first two lattice path- integrals. This will help clarify their physical meaning. 
For the two-step closed paths, and starting from any initial site, the electron retraces its first step on one of the six 
bonds connecting the initial site with its adjacent sites. This process can be designated symbolically by (•«-►), where 
the dot (•) indicates the initial site. The flux enclosed, of course, is zero. From these six closed paths of two steps 
each, we obtain 

S 2 = 6 (•<->•) = 6e 4 ° = 6 = z, 

where z stands for the coordination number of the cubic lattice. 

For the four-step closed paths, we need to consider four different possibilities: (1) the electron retraces twice each 

one of the six bonds connecting the initial site with its nearest-neighboring sites (•«-►); (2) after retracing once one of 
the six bonds connecting to the initial site, the electron retraces two steps on one of the other five bonds (<-*■•<-*■); (3) 
hopping first to one of the six adjacent sites, the electron retraces once one of the other five bonds (the one connecting 

to the initial site is excluded) and then returns to the original site (• < ~m > ); and (4) the electron traverses either 
counterclockwise or clockwise on one of the four elementary square cells (connecting to the initial site) on the yz-, 

zx- and xy-planes respectively [(•□ + •D) yz + (-n + •□) 22 , + (•□ + • □) a . 2/ ]. Thus 



S 4 = 6(-^)+30(< 
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(•□ +'°)yz + (•>=>+ + + •□)*„ 

= (6 + 30 + 30) e m 

+ 4 [(e w + e~ ia ) + (e 4b + e~ ib ) + {e tc + e~ ic )] 
= 66 + 8 (cos a + cos b + cos c) . 

To compute S21 in a systematic manner, we first define the quantity Sp} q , r , which is the sum over all possible paths 
of t steps on which an electron may hop from the origin (0,0,0) to site (p,q,r). From the definition of Sp} q ^ r , it is 
straightforward to construct the following recurrence relation for 



c(*+i) 
p>q? r 



exp ±i 



qc-rb \ s(t) 



p±l,q,r + eX P ±* 



ra — pc 



S, 



p,9±l,r + CX P( ^ 
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(3) 



This equation states that the site (p, q, r) can be reached by taking the (t+ l)th step from the six nearest-neighboring 
sites. The factors in front of the S^'s account for the presence of the magnetic field. By recursively using Eq. (3), 
we obtain the recurrence relation for Sp^ as 



S P,q? = 6 S p.q.r + exp[±i (qc - rb) } S^l 2 q r + cxp[±i (ra - pc) } S { £ q±2 r + exp[±i (pb - qa) ] S p t} q r±2 



+ 2 cos (|) [cx P ( ± ^ fl - b > 2 ^ 



+ 2 cos 



exp ±i 



. q(c — a) — (r — p)b 



S. 



(t) 
>± 

(«) 



.r(a + b) - (p + q)c\ (t) 



p±i, 9 ,r±i +expl =R 



q(c + a) - (r + p)b 



S. 



(t) 

p±l,(7,r=Fl 



+ 2 cos © [exp( ± ,fc^i^>) ^ ±1 , r±1+ ex P ( T , 



. p(b + c) - (q + r)a \ (t) 

2 ) S p,q±l.r T l 



(4) 
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We will now use Eq. (4) as the basic recurrence relation in our iteration scheme. Without loss of generality and 
for convenience, we choose the origin (0, 0, 0) to be our starting site. The initial conditions then read Sqq = 1 and 
Sp°l, r = for other (p, q, r)'s. It is evident that S21 is just equal to S^q . Putting t = in Eq. (4), we directly have 

"0,0,0 — °J 



0(2) 


c(2) 
— "1,-1,0 


o(2) 
— "-1,1,0 


c(2) 
— "-1,-1,0 


= 2 


o(2) 

1,0,1 


c(2) 
— "1,0,-1 


c(2) 
— "-1,0,1 


c(2) 
— "-1,0,-1 


= 2 


c(2) 
"0,1,1 


c(2) 
— "0,1,-1 


c(2) 
— "0,-1,1 


c(2) 
— "0,-1,-1 


= 2 



■(f) 



and 



q (2) _ q (2) q (2) _ q (2) _ q (2) _ q (2) _ -, 

"2,0,0 — "-2,0,0 — "0,2,0 ~ "0,-2,0 ~~ "0,0,2 ~~ "0,0,-2 — L - 



Using the above results for the Spar's and putting t = 2 in Eq. (4), we then obtain the S^/s. In general, by following 
this procedure we can obtain to any desired I. The properties discussed below will make the computation of 

Sq 2 q quite efficient. 

For a given I, a nonzero «Sp, 9 ,r exists only on those p, q, and r satisfying \p\ + \q\ + |r| = 0, 2, . . . , 21. This stems from 
the fact that only the sites (p, q, r) satisfying this condition can be reached by an electron after 21 hops from (0, 0, 0). 
The total number of these sites is 

(21) = 16Z 3 + 24Z 2 + 14/ + 3 
T 3 

Among the Sp i9 , r 's, we find the following symmetries hold for any I. We omit the superscript (2/) below. 

Sp,q,r{ a i b, c) = >S'_p._q._ r (a, 6, c), 



S—p,q,r [Oi, b, c) 


— Sp,-q. 


- r (a,b, 


c) 


*5p, — q.r {&, b, c) 


= S-p,q. 


- r (a,b, 


c) 


^p,q>— r bj c) 


— S— p t — 


q.r ( a , b : 


c) 



<Sp,<j,r(o, b, c) — »S— p,<j,r( a, 6, c) — "Sp, — 9,r( a , c ) — Sp, q _ — r (a,b, c), 

and 

Sp,q,r(.0'i b, c) = Sp,r,<j( a ; C, 6) = S ry p y q(b, C, d) = Sq : p t r{b, CI, c) = S q _ r ^p(c, CI, 6) = S r , q ,p(c, b, fl). 

We can therefore reduce our calculation to the S^^'s at sites (p,q,r) in the first octant with p > g > r > 0. Thus, 

(21) 

the total number of the independent 

(20 _ I 3 + 6l 2 + 12? + 9 - n 
1 ~ 9 

for I = 3m + n, where m are non-negative integers and n = 0, 1, 2. Thus the number of the S^'s to be computed is 
significantly reduced (about a factor of 48) compared to that of the whole Sp 2 q] r . 

Finally, it is worthwhile to note that to obtain the lattice path- integrals up to S2L, it is sufficient to compute the 
Sp^ys for < p + q + r < L when L is even and the S p f q ^'s ior0<p + q + r<L — 1 when L is odd. To be more 
specific, in doing the iteration of S^ q \ r for L/2 < I < L, we need only to concern ourselves with those S^r's with 
p + q + r = 0,2,...,2(L-0- 
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C. Results for lattice path-integrals with general flux orientations 

For the general flux orientation (a, b, c), we have obtained £2; up to 21 = 20. Here: J2( a ) denote sums over a — a,b, c; 
S( Q/ 3) denote sums over {a [3) = (ab), (be), (ca); and ^(a/37) denote sums over (a/37) = {a be), (bca), (cab). Also, 
for instance, the term cos(a ± (3) means cos(a + (3) + cos(a — (3). Below we present the results for S4, Sq, S$, and Sio- 

S 4 = 66 + 8 

cos a, 

(a) 

Sa = 876 + 240 cos a + 24 cos 2a] + 24 ^ cos(a ± (3) + 12 cos(a + b + c) + 12 ^ cos(a + /3- 7 ), 

(a) (a/3) (a/37) 

5 8 = 12978 + [ 5632 cos a + 1000 cos 2a + 96 cos 3a + 16 cos 4a ] 

(a) 

+ ^ [1120 cos(a ± 13) + 112 cos(2a ± /3) + 112 cos(a ± 2(3) + 32 cos(2a ±2/3)] + 576 cos(a + b + c) 

(a/3) 

+ [ 576 cos(a + /3 - 7) + 64 cos(2a + /3 ± 7) + 64cos(2a - /3 ± 7) 

(a/37) 

+ 16 cos(a + 2/3 ± 2 7 ) + 16 cos(a - 2/3 ± 2 7 ) ] , 
5i = 208836 + ^ [124080 cosa + 30040 cos 2a + 5040 cos3a + 1160 cos4a+ 160 cos5a + 40 cos 6a] 

(a) 

+ ^ [ 36680 cos(a ± (3) + 6800 cos(2a ± (3) + 6800 cos(a ± 2/3) + 2160 cos(2a ± 2/3) 

(a/3) 

+ 600 cos(3a ± (3) + 600 cos(a ± 3/3) + 240 cos(3a ± 2/3) + 240 cos(2a ± 3/3) + 80 cos(4a ± (3) 
+ 80 cos(a ± 4/3) + 60 cos(2a ± 4/3) + 60 cos(2a ± 4/3) + 40 cos(3a ± 3/3) ] + 19860 cos(a + b + c) 

+ [ 19860 cos(a + (3 - 7) + 4040 cos(2a + (3 ± 7) + 4040 cos(2a - /3 ± 7) 

(a/3 7 ) 

+ 1280 cos(a + 2(3 ± 2 7 ) + 1280 cos(a - 2(3 ± 2 7 ) + 400 cos(3a + (3 ± 7) + 400 cos(3a - (3 ± 7) 
+ 160 cos(a + 2(3 ± 37) + 160 cos(a - 2(3 ± 37) + 160 cos(a + 3/3 ± 2 7 ) + 160 cos(a - 3/3 ± 2 7 ) 
+ 40 cos(4a + (3 ± 7) + 40 cos(4a - (3 ± 7) + 40 cos(3a + 2(3 ± 2 7 ) + 40 cos(3a - 2(3 ± 2 7 ) 
+ 40 cos(a + 4/3 ± 2 7 ) + 40 cos(a - 4/3 ± 2 7 ) + 40 cos(a + 2(3 ± 4 7 ) + 40 cos(a - 2(3 ± 4 7 ) 
+ 20 cos(4a + 2/3 ± 2 7 ) + 20 cos(4a - 2(3 ± 2 7 ) + 20 cos(a + 3/3 ± 37) + 20 cos(a - 3/3 ± 37) ] . 

We have also worked out a considerable number of lattice path-integrals (up to the 40th order) for various special 
orientations of the flux. The results for S4, Sq,. . ., S12 are listed in Table I. The results for Su, S15,. . ., S40 are not 
presented. 



III. TOTAL KINETIC ENERGY 



We now proceed to calculate the total kinetic energy Et of the half-filled Fermi sea. We do this by using a series 
expansion of the total energy. We also present an alternative approach based on a moment expansion of the density 
of states. Through these two different approaches, we show that our goal here — the calculation of St(B) — is reduced 
to the evaluation of the lattice path- integrals £2; ■ 



A. Series expansion for the total energy 



Let us work on the {IV'i)} basis. At half-filling, the total kinetic energy of nonintcracting spinless electrons is the 
sum of the lowest N/2 eigenvalues, where TV is the total number of sites. Since the energy spectrum of the cubic 
lattice is symmetric under {E} — > {— E}, we can write the total energy per site as 



Ex 



E<0 



Ho 

z 



(5) 
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Here Tr denotes the trace and Hq is the corresponding diagonalized Hamiltonian of H. Notice that the absolute 
value is typically defined for scalar numbers. In this case, \Hq/z\ refers to an operator obtained by taking the absolute 
value of every matrix element of the operator Ho/z. Noting that —I < \Ho/z\ < I, we expand \Hq/z\ into a series in 
terms of Chebyshev polynomials T2k(H /z) as 



- I + - V W r T 2n ( 

7T 7T » 



ir ' An 2 - 1 

n=l 



where / represents the identity operator. Exploiting the equality 



V Z / ; _n 



(6) 



(7) 



where 



we obtain from Eqs. (5-7) 



(n + Z-1)! 
(2/)!(n-01' 



E 



N ^ 4n 2 

n=l 



X> (Trtf 2 ') 



(8) 



where we have replaced Tr(i7g') by Tr(H 21 ), as they are equal to each other. 
Assuming periodic boundary conditions on the lattices, we have 

TrH 2l =N{^\H 2l \^). 

The total kinetic energy per site is then given by 



(9) 



E T (a,b,c) = -- 

7T 



vZ=0 



(10) 



where we have assumed the highest order of the lattice path-integral obtained is 2L. The above result is exact when 
L = oo. Truncations at high values of L provide excellent approximations to the L = oo case, because the terms 
corresponding to L larger than either 4, 5, or 6 (depending on the orientation of the field) are negligible. Our results for 
the general field orientation (a, b, c) of Et include all terms up to L — 10. For the special flux configurations (0, 0, <j>), 
(0, (j>,4>), {(j),<j),4>), (0,0,-20), and (n,ir,<f>), we obtain all terms up to L = 20. Terms up to L = 10 (corresponding 
to 520 ) include contributions originating from ~ 10 14 three-dimensional closed paths, while terms up to L = 20 
(corresponding to S40) include contributions coming from ~ 10 29 closed paths in a 3D cubic lattice. Each one of these 
closed paths is weighted by its corresponding magnetic phase factor. 

Note that, by construction, the lattice path-integrals are local quantities and are valid for any value of <p. The 
periodic boundary conditions make the lattice path-integrals homogeneous (i.e., translationally invariant) and are 
unrelated to the imposition of a Bloch theorem. In contrast to most works studying tight-binding electrons in a 
magnetic field, here we never invoke a k— space or reciprocal space: our sums are all defined in direct space. 



B. Moment expansion for the density of states 



Here we analytically apply the method of moment expansion to obtain the total energy Et- The starting point 
now is the one-particle Green's function Ga(E) = (il>i\(E — -H") -1 ^), which can be expressed as 



00 

Gu(E) = - £ 



E 2l 



w OO 

= -Y 

E f-^ 



1=0 1=0 
The nonzero 2Z-th moment M21 of the density of states p(E) is given by 



52; 
W 



(11) 
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M 2 l = 



E 



P {E) dE. 



(12) 



Taking into account the relation between p(E) and the imaginary part of the Green's function, it can be derived that 

<?2Z 



After the change of variables u) = E/z, Eq. (12) can be rewritten as 

M 2 i = J u 21 p(u) duj, 



(13) 



(14) 



where p(u) — z p(ujz). 



We now expand p(uj) in a scries of Chebyshev polynomials weig hted by l/s/l — uu 2 as 

T 2 fc(a>) 



tr n VI- 



k=0 

Substituting Eq. (15) into Eq. (14) and using for I > k 

r 1 L0 2l T 2k (uj) 



1 Vi^P 



■ dui = L 



2l.2k 



(2l)\n 



4' (i - k)\ (i + ky: 



we obtain 



M21 = E C 2 fe hi,2k ■ 



(15) 



(16) 



(17) 



k=0 



Given the moments up to M2L, the coefficients C 2 k, for k = 0, . . . , L, can then be exactly determined one by one in 
terms of the moments. For instance, 



C = 



Mq = 1 

^0,0 71" ' 



and 



C 4 = 



I2.2 37r' 
M.4 — Co ^4,0 — C2 I42 



14,4 



4 
27^ 



1 + - (cos a + cos b + cos c) 



In general, every C2 k can be computed from the previous ones by 

1 / k ~ 1 

C 2 fe = 7 M-2k — > C21 hk,2 

hk,2k \ ~ 



1=0 



Noting that at half-filling the total energy Et is 
and using the result 



( Ep{E)dE, 

J —z 



(18) 



(19) 



u)T 2 k{u) 



dui = 



(-!)< 



yr^2 (2k-i)(2k + i) : 
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wc therefore obtain the following expression for the total kinetic energy 



L 

E T (a,b,c) = z ^2 

fc=0 



(2k- l)(2Jfe + l) 



C 2 k(a,b,c). 



(20) 



It is interesting to notice that Eqs. (10) and (20) lead to the same result for the total energy Ex, even though we 
used quite different approaches. Both results are exact when L = oo. However, terms with k larger than 4, 5, or 
6 (depending on the field orientation) provide negligible contributions to the sum. Our results for the general field 
orientation (a, b, c) of Et(cl, b, c) include all terms up to L = 10, while for the five special flux configurations studied 
below (Sec. III.C), we include all terms up to L — 20. 



Substituting the expressions for the 6>2z(a, b, c)'s into Eq. (10) or the results for the C2/(a, b, c)'s into Eq. (20), we 
find that the lowest energy state is reached when a = b = c = tt (namely, a half flux quantum per plaquette on the 
yz-, zx- and xy-planes) . We thus provide an explicit analytic derivation for the absolute minimum energy of the flux 
state. The result is consistent with a theorem proved by Lieb.cl 

We also obtain the total energy Et(4>) as a continuous function of the flux for the following flux orientations: 
(0,0,0), (0,0,0), (0,0,0), (0,0,-20), and the asymmetric case (tt,tt, 0). Notice that the field direction of (0,0,0) 
is perpendicular to that of (0, 0, —20). These results are plotted in Fig. 1. They are calculated through Eq. (20) by 
using their respective CWs for 2k = 0, 2, 4, . . . , 40. Thus, we have added the contributions of ~ 10 29 electron closed 
paths in a 3D cubic lattice, each one weighted by its magnetic phase factor. 

For comparison purposes, we also present E^ D \(j)) for the 2D case, obtained by using the corresponding C^t s 

up to C^ D \ This corresponds to summing over ~ 10 44 electron closed paths in a 2D square lattice, each one with its 
corresponding phase factor. Note that the absolute minimum of Et(4>) occurs at = tt in all of these flux orientations, 
except for the orientation (0, 0, 0) in 3D. It becomes clear that hopping in an extra dimension drastically changes the 
properties found in strictly two-dimensional systems. 

In Fig. 2, we plot the total kinetic energy of electrons at half- filling Et for the 3D cubic lattice for various field 
orientations versus the order of the highest-order lattice path-integral used to calculate them. For comparison, we 
also show the numerical values of Et obtained in Ref. 5, which mostly has only three significant figures. The Et 
reach steady values which are consistent, within ~ 2%, with numerical ones for lattice path-integrals of order 2L equal 
to either 4, 6, or 8, depending on the orientation of the field. Higher order lattice path-integrals provide negligible 
contributions to Et- For instance, for (0,0,0) both approaches give the same result, within ~ 0.5%, for 2L > 4. 



The approach described in section III.B can be directly generalized to calculate the total energy Et(v) for any 
electron filling z;(0<^<l/2)as 



C. The lowest-energy flux state and Et for various flux orientations 



D. Total energy for any filling 




(21) 



where the Fermi energy \i is determined by 




(22) 



Utilizing Eq. (15) and carrying out the integrals, we obtain 




(23) 



where 9 = arccos(^/z) (tt > 9 > n/2), and given u, 9 can be solved from the following equation 
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V = 1 



fe=i 



2A-" 



sin(2fc(9) 
~2fc ' 



(24) 



Note that < v < 1/2 corresponds to -z < \i < and tt > 9 > n/2. Writing v = 1/2 - 8 and 6> = tt/2 + 77, Eqs. (23) 
and (24) can then be rewritten as 



E T ( - - 8 



z 

IT 



cosr7 + !^(-l) fe C7 2fc 



fc=i 



i[{2k - 1)77] cos[(2fc + 1)?/] 



2fc- 1 



2k + 1 



and 



fc=l 



sin(2fc?7) 
2fc ' 



(25) 



(26) 



In the small S limit (i.e., the electron filling v is close to one half and 77 is very small) sin(2fc?7) ~ 2fcfy; we then have 

7T 



V 



1 + 7rC(a, 6, c) 



5, 



(27) 



where 



C(a,6,c) = ^(-l)' £ C 2fe (a,6, C ). 



fe=i 



Therefore, using the approximation cos[(2fc ± 1)77] ~ 1 — (2k ± l) 2 ?7 2 /2, we obtain the total energy £r(l/2 — 8) close 
to (and below) half-filling 



E T ( - - 8 



7r V 1 + irC 



8)+z 



1 — cos 



fc=i 

7T 



V (-l) fc r * U 

^ (2/fc - l)(2fc + 1) 2fe 2 (1 + to 



(1 + ttC) 5 



1 + ttC 



2(1 + ttC) 2 



(28) 



Finally, by differentiating both Et in Eq. (23) and v in Eq. (24) with respect to /i, we can establish the following 
identity 



A,Et dv 



(29) 



It is worthwhile to emphasize that the formalism described in this section is equally applicable to the 2D square lattice 
case. 



IV. MAGNETIC MOMENT AND ORBITAL SUSCEPTIBILITY 



In this section we investigate two experimentally accessible observables, namely, the magnetic moment M and the 
orbital susceptibility \ °f this quantum system. It is well knownD that at absolute zero temperature M and \ arc 
given by the first- and second-order derivatives of the total energy 

« = -2| (30) 

and 

X = -^&- (31) 

With our analytical results for Et at hand, equations (30) and (31) provide a straightforward way to study these 
important quantities for any orientation of the externally applied field. For illustration purposes, below we will focus on 
the following flux orientations: (0, 0, cf>), (0, </>, </>), (</>, <fi, cf>), and (</>, 4>, —2(f)). It will be seen that the magnetic response, 
in the presence of a strong periodic potential, is significantly different from the familiar Landau diamagnetism of a 2D 
electron gas, where — \ has its largest value at B = and its vicinity and decreases monotonically as B is increased. 
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A. Zero-field susceptibility x(/ i ) = X{4> — 0, fJ-) versus Fermi energy /j, 



When we increase the electron filling factor is, the Fermi energy /x increases. We therefore first study the zero-field 
susceptibility x(m) = x(0 — 0,/i) as a function of fi through 



^From Eq. (23), we readily obtain 



xO) 



d 2 E T 



dB 2 



(32) 



B=0 



x(m) 




d 2 C 2k &) 



0=0/ 



sin[(2fe - 1)6»] sin[(2A: + 1) 



2k- 1 



2k + 1 



(33) 



where d 2 C2k{4>) / d4> 2 \ ( p=Q can be easily computed for each of the flux orientations under consideration. 

In Fig. 3 we plot the negative of the susceptibility — x as a function of the Fermi energy /x for various flux configu- 
rations. Large diamagnetic (i.e., negative) susceptibilities are observed at very low electron filling. For increasing /z, 
the quantity — x decreases and fluctuates around zero. For /i larger than a certain value (~ —3.2), paramagnetism 
prevails. Indeed, the orbital response is weakly paramagnetic at half-filling (fi — 0). These results are in qualitative 
agreement with the observations of Skudlarski and VignaleQ based on calculating the current-current correlation func- 
tion. Furthermore, note that the absolute value of x f° r (0, 0, 0) ls always twice that for (0,0, 0); |x| for (0,0,0) is 
always three times that for (0, 0, 0); and |x| for (0, 4>, — 20) is twice that for (0, 0, cf>). These results clearly reflect the 
orbital origin of the susceptibility. 



B. Field-dependent M(0) and x{4>) at v = § 

We now focus on the system at half-filling and examine the continuous dependence of the magnetic moment M and 
susceptibility x 011 the magnetic flux 0. Employing Eq. (20) for £t(0), we can directly obtain 

mi dE T (4>) A (-l) k dC 2fc (0) 
— M{<fj) = — n -= z ^ (2fc _ ^ (2fc+ i} (34) 



/c=0 



and 



d 2 E T ^) ^ (-l) fc d 2 C 2k (cb) 
X{(t>) = = Z ^ (2A-l)(2fe+l) (35) 

The results for four different flux orientations are shown in Fig. 4. Let us look at —M{<jj) — —M{(j),v — 1/2) first. 
Starting from = 0, —M{<p) decreases from zero and remains negative with some fluctuations. There are then 
some irregular "sine-function-shaped" oscillations present. These low-field oscillations in —M{4>) are observed for 
all orientations shown in Fig. 4. For the (0, 0, (j>) flux orientation, — M has a positive value for (f)/2n k, 0.3-0.5, 
reaching zero at <j) = tt. On the other hand, for (0, 4>,<f>), (0, </>,</>) and (0,0,-20), —M is mostly negative for low 
fields (0/27T « 0-0.1), oscillates around zero for intermediate fields (0/27T « 0.1-0.3), is negative for larger fields 
(0/27T w 0.3-0.5), and reaches zero at = it. At relatively large flux values (0/27T « 0.4-0.5), — M: decreases to 
zero (at = tt) from a positive value for (0, 0, 0); and increases to zero (at = 7r) from a negative value for (0, 0, 0), 
(0,0,0), and (0,0,-20). 

The field-dependent orbital magnetic susceptibility, — x = — x(0: ^ — 1/2), fluctuates somewhat evenly around 
zero in the low-field regime. However, — x fluctuates less, and around a small negative value for the flux orientation 
(0,0,0), where only one plane (xy-plane) is penetrated by the flux. The fluctuations are more pronounced in the 
flux orientations: (0,0,0), (0,0,0), and (0,0,-20), where at least two perpendicular planes are affected by the 
field. As the flux is raised, — x tends to fluctuate less and: (i) a crossover from diamagnetism to paramagnetism is 
observed in the configuration (0, 0, 0); and, on the other hand, (ii) a crossover from paramagnetism to diamagnetism 
is observed in the orientations (0, 0, 0), (0, 0, 0), and (0, 0, —20). In other words, at = tt, the magnetic response is 
weakly paramagnetic in the flux configuration (0, 0, 0), while the other three orientations provide a strong diamagnetic 
response in the lattice. 

Figure 4 suggests that, for both M and Xi the following four quantities tend to grow larger with the increase of 
planes exposed to the field: the number of oscillations, the frequency of the oscillations as a function of the flux, the 



10 



number of nodes, and the amplitude of the oscillations. For instance, M(<j>, </>, <ft) oscillates more rapidly and strongly 
than M(0, <f>, 4>), which at the same time oscillates more rapidly and strongly than M(0, 0, <fi). A similar general trend 
is observed between x{4>i 4>: 4>), x(0, 0, 0), and x(0,0, </>). These features can be understood intuitively: the more 
perpendicular planes are exposed to the flux, the more strongly M and \ wm be affected. 

V. SUMMARY 

In conclusion, we present an investigation of quantum interference phenomena of tight-binding electrons on the 
3D cubic lattice in a continuously-tunable magnetic field with arbitrary orientation. Previous work on this problem 
focused on a discrete set of (rational) magnetic field values. We study the total kinetic energy, and subsequently the 
magnetic moment and orbital susceptibility. The main results include: (1) an analytic study of electron quantum 
interference effects resulting from sums over magnetic phase factors associated with 3D closed paths, (2) a very 
efficient computation of these "lattice path-integrals" S21 in closed-form expressions, (3) explicit analytic expressions, 
in terms of the lattice path- integrals, for the Fermi-sea ground-state energy Et as a function of the fluxes for electron 
fillings v at and near one half, (4) the 3D lattice path-integrals to very high order and the total energies in various 
flux orientations, (5) the investigation of the zero-field orbital susceptibility xip) as a function of the Fermi energy /z, 
and (6) the magnetic moment M((f>) and susceptibility x(</>) as functions of the flux, both at half-filling. 

We find that the absolute minimum of Et(4>) at half-filling occurs at <p = it in all of the flux orientations under 
consideration, except for the configuration (0,0, (f>) in 3D. It becomes evident that hopping in an additional direction 
drastically changes the properties found in strictly two-dimensional systems. It is also seen that the magnetic response, 
in the presence of a strong periodic potential, is significantly distinct from the familiar Landau diamagnetism — where 
in a 2D electron gas —\ takes the largest value close to B = 0, and decreases monotonically with increasing B. 
For the zero-field susceptibility x{^ JL )^ diamagnetism dominates in spite of small fluctuations of xif) around zero for 
small electron filling v. On the average, with increasing /i, the quantity —xip) non-monotonically decreases from 
a large positive (diamagnetic) value to a relatively small negative (paramagnetic) value. For /i > —3.2, the orbital 
response becomes paramagnetic. Both the field-dependent —M(<f) and — xi'P) exhibit irregular oscillations according 
to the direction of the field. For the four flux orientations (0, 0, 0), (0, <j>, (j)), (0, <j>, </>), and (0, <j>, —2(f)), the magnetic 
moment M((f>) is always zero at 4> = and tt; and paramagnetism (\ > 0) exists at cj> = for all these flux orientations. 
However, when 0/27T = 1/2, the magnetic response is paramagnetic for the flux configuration (0, 0, (/>), and diamagnetic 
for the other three orientations. _ 

Here, we add a remark on the quantity dET(4>)/d<fi. It is known that at zero temperature the persistent currents 
in a metal ring threaded by a magnetic flux tp is proportional to the sum over the contributions of dE n /dip from all 
occupied states, where E n is the eigenstate energy. We therefore can regard dET{<f>) / d(f> as a generalized "current" in 
this multiply-connected lattice system. 

An approach commonly employed in recent years to study electrons in a magnetic field <p = p/q (e.g., Refs. 4-6) 
uses the Bloch theorem and maps the problem into a q x q matrix problem — related to a q x q cell with periodic 
boundary conditions (PBC) in the actual lattice. Thus, for <f> — 1/2, the electron energy levels are determined by 
considering a 2 x 2 cell with PBC in the lattice and diagonalizing a 2 x 2 matrix. This approach presents a problem 
for any irrational field, since q — > oo, and a periodic cell cannot be realized. We do not follow this approach, thus, we 
can explicitly consider a continuum of values of 4>, without treating irrational numbers on a special footing. 

Our lattice path- integrals are local quantities. By construction, they are valid for any value of the magnetic 
field. The PBC we impose [in Eq. (9)] are only designed to make these lattice path-integrals homogeneous (i.e., 
translationally invariant) and have nothing to do with the imposition of a Bloch theorem. Indeed, in contrast to most 
works studying electrons in a magnetic field, here we never invoke any k-space or momentum-space: our sums over 
closed paths (lattice path-integrals) are all explicitly defined in real space. 

The lattice path-integrals obtained here are "many-loop" generalizations of the standard "one-loop" Aharonov- 
Bohm-type argument, where the electron wave function picks up a phase factor e 1 * each time it travels around a 
closed loop enclosing a net flux <I>. The evaluation of these lattice path-integrals enables us to analytically obtain 
the total energies, magnetic moments, and orbital susceptibilities of the corresponding flux states. The spirit of our 
approach follows Feynman's programme: to derive physical quantities in terms of "sums over paths" . This method 
is considerably different from the standard ones that have been employed so far (e.g., transforming the problem 
to momentum space and computing Et numerically). In particular, it allows the analytic calculation of physical 
quantities as explicit functions of a continuously-tunable flux, while other approaches need to separately consider 
different cases (e.g., matrices) for several discrete (rational) values of the magnetic field. The lattice path approach 
can also be used for a variety of other physical problems, including the derivation and analysis of the superconducting 
transition temperature in wire networks and Josephson-j unction arrays (see, e.g., Refs. 10 and 11), and the analytical 
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computation of the magnetoconductance for strongly localized electrons. 
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TABLE I. Lattice path-integrals S21 (with the order 21 — 4, 6, 8, 10, 12) for the 3D cubic lattice in various flux configurations: 
(0, 0, 0), (0, 0, 0), (0, 0, 0), (0, 0, —20), and the asymmetric case (tt, it, 0). 

21 S 2 i in (0,0,0) 

~4 82 + 8 cos 

6 1452 + 384 cos + 24 cos 20 

8 29218 + 13440 cos + 1960 cos 20 + 96 cos 30 + 16 cos 40 

10 638756 + 422800 cos + 96760 cos 20 + 1 1760 cos 30 + 2284 cos 40 + 160 cos 50 + 40 cos 60 

14865220 + 12737640 cos + 3909984 cos 20 + 764216 cos 30 + 182328 cos 40 
+28392 cos 50 + 7680 cos 60 + 528 cos 70 + 144 cos 80 + 24 cos 90 

21 S-2i in (0,0,0) ~ 

4 74 + 16 cos 

6 1188 + 576 cos + 96 cos 20 

8 22186 + 16736 cos + 4976 cos 20 + 736 cos 30 + 96 cos 40 

10 460236 + 461760 cos + 186360 cos 20 + 49760 cos 30 + 11840 cos 40 + 2080 cos 50 + 520 cos 60 
10354220 + 12623616 cos + 6211104 cos 20 + 2269840 cos 30 + 742800 cos 40 

+209328 cos 50 + 67008 cos 60 + 14016 cos 70 + 3744 cos 80 + 480 cos 90 

21 S21 in (0,0,0) 

4 66 + 24cos0 

6 948 + 756 cos + 144 cos 20 + 12 cos 30 

8 16626 + 19392 cos + 6744 cos 20 + 1584 cos 30 + 336 cos 40 + 48 cos 50 

10 338616 + 483420 cos + 230340 cos 20 + 82980 cos 30 + 27000 cos 40 + 7740 cos 50 + 1980 cos 60 + 420 cos 70 + 60 cos 80 

7672212 + 12250440 cos + 7095780 cos 20 + 3290664 cos 30 + 1370952 cos 40 + 528120 cos 50 + 193524 cos 60 

+65736 cos 70 + 20952 cos 80 + 5952 cos 90 + 1512 cos 100 + 288 cos 110 + 24 cos 120 

21 S21 in (0,0,-20) ~ 

4 66 + 16cos0 + 8cos20 

6 912 + 528 cos + 336 cos 20 + 48 cos 30 + 36 cos 40 

8 14930 + 14016 cos + 10080 cos 20 + 3040 cos 30 + 1992 cos 40 + 384 cos 50 + 240 cos 60 + 32 cos 70 + 16 cos 80 

282796 + 356640 cos + 275560 cos 20 + 122640 cos 30 + 79780 cos 40 + 28160 cos 50 + 17280 cos 60 
+5280 cos 70 + 3040 cos 80 + 800 cos 90 + 440 cos 100 + 80 cos 1 10 + 60 cos 120 

6045696 + 9120912 cos + 7393704 cos 20 + 4172032 cos 30 + 2809680 cos 40 + 1342512 cos 50 
12 +841808 cos 60 + 367440 cos 70 + 220836 cos 80 + 90384 cos 90 + 52800 cos 100 + 19680 cos 1 10 

+11856 cos 120 + 3696 cos 130 + 2112 cos 140 + 576 cos 150 + 360 cos 160 + 48 cos 170 + 24 cos 180 

21 S2I in (7T,7T,0) 

1 50 + 8 cos 

6 492 + 192 cos + 24 cos 20 

8 5410 + 3456 cos + 808 cos 20 + 96 cos 30 + 16 cos 40 

10 64676 + 56720 cos + 18680 cos 20 + 4080 cos 30 + 1000 cos 40 + 160 cos 50 + 40 cos 60 

826820 + 900840 cos + 372384 cos 20 + 113336 cos 30 + 35448 cos 40 
+9192 cos 50 + 2880 cos 60 + 528 cos 70 + 144 cos 80 + 24 cos 90 
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FIGURE CAPTIONS 



FIG. 1. Lower frame: Total kinetic energy of electrons at half-filling Et{4>) for the 3D cubic lattice for various flux orien- 
tations: (0,0,0) (long dash), (0,0,0) (short dash), (0,0,0) (solid), (0,0,-20) (dot), and (n, n, 0) (dot-dash). Upper frame: 
E { T D \<t>) for the 2D square lattice in a perpendicular field. The flux values where the minima occur are indicated. We obtain 
these results by using the lattice path-integrals up to the order 40 for the 3D case and 76 for the 2D lattice — which corresponds 
to summing contributions over ~ 10 29 electron closed paths in a 3D cubic lattice and ~ 10 44 electron closed paths in a 2D 
square lattice. Notice that Et(4>) has the absolute minimum at 0/27T = 1/2 for all of these flux orientations, except for the 
configuration (0,0,0). 

FIG. 2. Total kinetic energy of electrons at half-filling Et for the 3D cubic lattice for various field orientations: (1/2, 1/2, 1/2), 
(1/3,1/3,1/3), (1/4,1/4,1/4), (0,1/2,1/2), (0,1/3,1/3), (0,1/4,1/4), (0,0,1/2), (0,0,1/3), and (0,0,0). Here (a,fo,c) means 
that the three fluxes through the elementary plaquettes on the yz-, zx- and xy-planes are respectively a/2n, b/2n and c/2n. 
The horizontal axis (2L) denotes the order of the highest-order lattice path- integral used in computing Et- The dashed lines 
indicate the corresponding values from Ref. 5, where only the first three significant figures are available, except for the cases 
(1/2, 1/2, 1/2) and (0,0,0). The Et reach steady values which are consistent, within ~ 2%, with numerical ones for lattice 
path-integrals of order 2L equal to either 4, 6, or 8, depending on the field orientation. Higher order lattice path-integrals 
provide negligible contributions to Et- Notice that the vertical axes here cover very narrow range of values (~ 10~ 2 ). For 
(1/2, 1/2, 1/2), (0, 1/2, 1/2) and (0, 0, 1/2), both approaches give the same result within ~ 1% for 2L > 8 and within ~ 2% for 
21 > 8, 8, and 4, respectively. For (1/3, 1/3, 1/3), (0, 1/3, 1/3) and (0,0, 1/3), both approaches provide the same result within 
~ 1% for 2L > 12, 10, and 12, and within ~ 2% for 21 > 8, 4, and 4, respectively. For (1/4, 1/4, 1/4), (0, 1/4, 1/4) and (0, 0, 0), 
both approaches yield the same result within ~ 1% for 2L > 16, 4, 4 and within ~ 2% for 21 > 4, 4, and 4, respectively. 

FIG. 3. The negative of the zero-field susceptibility — x(aO as a function of the Fermi energy \i for various flux configurations: 
(0,0,0), (0,0,0), (0,0,0), and (0,0,-20) corresponding to curves a, b, c, and d, respectively. Curve a is always the closest 
to the dotted reference line x — 0, while curve d is always the farthest from it. For all these orientations, the susceptibility 
exhibits a very non-monotonic behavior. For small electron filling v (i.e., small n), diamagnetism (x < 0) dominates in spite of 
small fluctuations of — x around zero. With increasing n, on the average, — x decreases from a large positive value to a negative 
one. For fi > —3.2 the orbital response is paramagnetic. 

FIG. 4. The negatives of the magnetic moment — M(0) (left column) and the orbital susceptibility — x(0) (right column) 
at half-filling as a function of the flux for various flux orientations as indicated. /.From top to bottom they are: (0,0,0), 
(0, 0, 0), (0, 0, 0), and (0, 0, —20). Both — M and — x exhibit non-periodic oscillations. Notice the change of the amplitude and 
the frequency of the oscillations as a function of the flux for different orientations of the field. At = and tv, the magnetic 
moments are always zero. Moreover, paramagnetism exists at = for all these flux orientations. However, when 0/27T = 1/2, 
the magnetic response is paramagnetic for flux configuration (0,0,0), and diamagnetic for the other three orientations. 
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